Nonuniversal intensity correlations in 2D Anderson localizing random medium 
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Complex dielectric media often appear opaque because light traveling through them is scattered 
multiple times. Although the light scattering is a random process, different paths through the 
medium can be correlated encoding information about the medium. Here, we present spectroscopic 
measurements of nonuniversal intensity correlations that emerge when embedding quantum emitters 
inside a disordered photonic crystal that is found to Anderson-localize light. The emitters probe 
in-situ the microscopic details of the medium, and imprint such near-field properties onto the far- 
field correlations. Our findings provide new ways of enhancing light-matter interaction for quantum 
electrodynamics and energy harvesting, and may find applications in subwavelength diffuse-wave 
spectroscopy for biophotonics. 

PACS numbers: (42.25.Dd, 42.25.Fx, 46.65.+g, 42.70.Qs) 



Correlations are of paramount importance in science and 
technology and offer fundamental insight into disparate dis- 
ciplines ranging from establishing the role of inflation in the 
early universe^, through mesocopic transport of electrons and 
photons in complex medicP, to the non-locality of quantum 
mechanics observed at the microscopic scaled In photonics, 
recording correlation functions provides a powerful method 
of revealing information of the quantum state of lig hP and 
correlations may be exploited for imaging beyond the clas- 
sical diffraction limhP. The emerging field of nanophotonics 
promises a number of new applications for controlling and 
manipulating light. However, unavoidable inhomogeneities 
lead to random multiple scattering that requires a statisti- 
cal description. Surprisingly, such multiple scattering can give 
rise to pronounced correlations between different propagation 
paths through the complex medium, at variance with the in- 
tuitive perception that random scattering fully scrambles all 
information. Such correlations can be either of classical or 
quantunP character and are traditionally universal, i.e., they 
depend solely on a single parameter, the universal conduc- 
tance, and are independent of the microscopic details of the 
medium. Placing a light source inside the medium creates a 
fundamentally new situation. In this case, nonuniversal corre- 
lations in the far field have been predicted^ that depend sen- 
sitively on the local environment of the embedded emitter^, 
and thus can be employed for probing the local properties of 
a complex medium. Consequently, recording the nonuniversal 
correlation function has been proposed as a novel method for 
in-situ spectroscopy with ultra-high resolution in a complex 
and disordered medium^, and algorithms have been proposed 
for imaging a source in an opaque medium from measurements 
of the local electromagnetic field density^. In this Letter, we 
measure for the first time the nonuniversal correlations gener- 
ated by single dipole emitters in a strongly scattering random 
medium. 

Figure [lja) illustrates an intensity speckle pattern gener- 
ated by a dipole source embedded in a disordered photonic 
crystal membrane where random noise in the position of the 
air holes gives rise to strong multiple scattering. The gener- 
ated correlations are gauged by the intensity correlation func- 
tion (denoted C o) that connects the intensities at two different 
spatial position j 8 * 9 * 12 ! ri and r 2 , 



where the brackets (. . .} nf y ff denote spatial averaging in the 
near- and far-field, respectively. Eq. |l]) holds after spectral 
averaging where all contributions from standard (i.e. univer- 
sal) mesoscopic correlations^ vanish, which is the case in the 
present experiment. The Co correlation has infinite range, i.e., 
it is independent of the observation points ri and r 2 in the far 
field. Eq. |l]) demonstrates that Co is linked to local near-field 
properties of the random medium through the local density 
of optical states p(ro) at the position ro. This provides a way 
of extracting Co for light since in optics, as opposed to, e.g., 
electronics, excellent point sources exist such as single atoms 
or quantum dots (QDs), and their spontaneous-emission dy- 
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Figure 1: (Color online) Intensity fluctuations in disor- 
dered photonic crystals, (a), Illustration of a far-field intensity 
speckle pattern originating from a single emitter in a mem- 
brane perforated by a disordered arrangement of holes. The 
red arrows illustrate the link between the near-field environ- 
ment of the emitter, the local density of optical states p(i"o), 
and the far-field intensity pattern, (b), (d), Scanning elec- 
tron micrographs (top view) of photonic crystal membranes 
with introduced disorder on the hole positions of S = 0% and 
S — 12%, respectively, (c), (e), QD photoluminescence spectra 
collected on the two samples. The spectra are scaled equally 
for better comparison. The bright peaks in the spectra are 
signatures of Anderson-localized modes. 
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Figure 2: (Color online) Spatial and spectral mapping of the far-field intensity pattern, (a), Scanning electron micrographs (top 
view) of photonic crystal membranes with a = 280 nm, r = 0.274a, and with the degree of disorder varying between 8 — 0% and 
8 = 12%. (b), Normalized high-power photoluminescence spectra collected while scanning the samples with a wavelength bin 
size of 1 nm. (c), High-power photoluminescence spectra collected while scanning photonic crystals with 8 = 6%, a = 280 nm, 
and different hole radii. The horizontal lines represent the spectral positions of the high-energy band edge of the unperturbed 
structures, (d), Normalized intensity distribution in a photonic crystal with 8 — 0% and 8 — 6%. The large deviation from 
Rayleigh statistics (dashed line) in the disordered case is due to Anderson localization. The solid red line represents the best 
fit to the theory of Ref. [6]. (e), Dimensionless conductance, g, and variance of the intensity fluctuations for various amounts 
of disorder. The lines are guides to the eye and the shaded area represents the samples where fully developed speckles were 
not observed for the sample size of 7 nm x 7 u.m employed in the present experiment. Anderson localization is observed for 
var(7) > 7/3 (dotted line). 



namics probes p(ro). Co correlations are distinct from other 
multiple scattering correlations due to the dependence on the 
local properties and this nonuniversality has been proposed as 
a method for subwavelength imaging of a random medium^. 

Probing the spatial fluctuations of p(i"o) with single emit- 
ters constitutes an efficient way of extracting Co. Founding 
experiments on modified emission dynamics in random me- 
dia employed ensembles of emitter j 14 " 16 ^, which do not en- 
able extracting p(i"o) in a complex medium and therefore 
the Co correlation could not be obtained. Experiments on 
single emitters have so far concentrated on cavity quantum- 
electrodynamics studies in one-dimensional disordered pho- 
tonic crystal waveguides^, where the collection of statistics 
required for obtaining the Co correlation could not be ob- 
tained. Here, we present measurements of the Co correlation 
in disordered photonic crystal membranes in the regime of 
Anderson localization^. We employ single InAs QDs that 



have recently been shown to be reliable probes of p(ro) in di- 
electric nanostructures when the internal fine-structure split- 
ting and the non-radiative processes of the QD are properly 
accounted for^H We observe the transition to the regime of 
two-dimensional Anderson localization by varying the degree 
of disorder of our samples, which gives rise to a pronounced 
Co correlation. Deep in the localization regime the Co corre- 
lation is found to be very strong, reaching a value as high as 
3.5, which is several times larger than is expected for scalar 
waves i n the m ost strongly scattering 3D structures reported 
to dat c 18 * 20 ^ . Finally, we demonstrate that the statistical 
properties of the Anderson-localized modes, and therefore the 
Co correlation, can be controlled and fine tuned by employing 
the underlying dispersion of the photonic crystal. The aid of 
a photonic band edge for Anderson localization of light was 
the original motivation of S. John for proposing photonic band 
gap crystals^, and our experiment provides to our knowledge 



3 




Figure 3: (Color online) Exploring the near-field by probing dynamic properties of single emitters, (a), Decay curves for two 
QDs at different positions in a photonic crystal with 8 = 6% and the fits to bi-exponential models (red lines). QD X is enhanced 
by the coupling to an Anderson-localized mode while QD 2 is inhibited, (b), Radiative decay rate (7 ra d) distributions measured 
at A = 910±3nm in the bulk membrane and in photonic crystals with a = 280 nm and r = 0.274a affected by different amounts 
of disorder. The yellow bars show the 7 ra d distribution measured in a photonic crystal with the same lattice constant but with 
r = 0.245a. (c), Co correlation function as extracted from the 7 ra< j distributions (black symbols) and from the spatial intensity 
fluctuations (red symbols). The shaded area represents the region of not fully developed multiple scattering and the lines are 
guides to the eye. (d), Variation of Co x ^fg with degree of disorder illustrating the nonuniversality of Co- The curve is a guide 
to the eye. 



the first direct experimental demonstration of this effect. 

To probe the far-field properties of our samples, we ex- 
cite many QDs embedded in the system and measure the 
far-field emission pattern. Fig. |TJb) and [TJd) show top-view 
scanning electron micrographs of photonic crystals where the 
holes have been randomly displaced corresponding to 8 = 0% 
and 8 — 12%, where 8 is the standard deviation of the hole 
positions with respect to the lattice constant. Fig. [ljc) and 
[lje) show recorded photoluminescence intensity spectra of the 
light emitted from the QDs under strong pump power condi- 
tions where the photonic modes of the sample are revealed, 
as opposed to the single QD emission lines seen at low excita- 
tion powei^. A smooth emission spectrum is observed for no 
intentional disorder (8 = 0%), which resembles the inhomoge- 
neously broadened emission spectrum of the QDs. Increasing 
the degree of disorder leads to larger intensity fluctuations and 
for 8 > 5% pronounced resonances are observed. This is due 
to the imperfections in the hole position that lead to random 
multiple scattering eventually inducing Anderson localization 
if the disorder is sufficiently pronounced. Figure[2ja) and[5|b) 
show a range of photonic crystals with various degrees of dis- 
order and the far-field intensity collected while raster scanning 
the samples. By varying either the hole radius (Fig. I2fc)) or 
the lattice constant (not shown) for a fixed amount of dis- 
order (8 = 6%), we observe that the spectral region where 
the cavities appear follows the photonic crystal band edges, 
which in turn implies that the Co correlations can be tuned 
by taking advantage of the underlying dispersion, as will be 
discussed later. The Anderson-localized cavities appear in a 
spectral range close to the band edg«P (A = 903 nm). This 
range broadens with the amount of disorder, which is a signa- 



ture of the so-called Lifshitz tail that was originally proposed 
for electrons^ and originates from the broadening of the band 
edge with disordei^. 

The signatures of the nonuniversal Co can be demonstrated 
by extracting from the far-field intensity measurements the 
universal dimensionless conductance g. g is the universal scal- 
ing parameter that governs the statistical aspects of light 
transport in a random mediurrP^, and the threshold for local- 
ization in absence of loss is^ g = 1 . g can be obtained from the 
probability distribution describing the intensity fluctuations® 
P(I = 1/(1)), which is extracted experimentally from the in- 
tensity scans shown in Fig. [2jb), and displayed in Fig. ^d). 
See Supplementary Information for further details. In the 
highly disordered case, the distribution significantly deviates 
from the Rayleigh statistics expected for diffusive propaga- 
tion. Excellent agreement with the theory of Ref. 6 is ob- 
served over five orders of magnitude, which is quite remark- 
able since the theory is formally derived only in the meso- 
scopic regime at the onset of Anderson localization. Here we 
find it valid also deeply inside the localized regime in agree- 
ment with previous experimental work on ultrasouncP^. From 
the comparison with our data we obtain g = 0.22 ± 0.07 for 
8 = 6%. We note that pronounced loss may reduce the value 
of g, but cannot explain a large variance. In the presence of 
loss, var(J) > 7/3 serves as a localization criterion for a mul- 
timode waveguide^ that can be extracted directly from the 
experimental data without assuming any model. Note that for 
this criterion it is assumed that the intensity fluctuations of 
the light leaking vertically out of the structure obey the same 
statistics as the transmission. Numerical modeling of ID sys- 
tems has proven the validity of such an assumption, and the 
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required variance of intensity fluctuations for Anderson local- 
ization was in fact found to be smaller than the transmission 
fluctuations^. The Anderson-localization criterion is clearly 
fulfilled in our experiments, and we observe, e.g., var(J) = 7.2 
for S = 6%. Figurej^e) displays g and var(J) versus the degree 
of disorder and shows that an optimum degree of disorder ex- 
ists (S = 6%) where g is smallest and var(7) largest, thus light 
localization strongest. Such an optimum has been predicted 
theoretically for 3D disordered photonic crystals^, and re- 
flects the fact that disorder is introduced as a perturbation 
to an underlying ordered structure that aids Anderson local- 
ization. For large degrees of disorder, the 2D band gap closes 
and band edge localization disappears. 

The near-field environment is investigated by probing dy- 
namic properties of the embedded emitters. For weak light- 
matter interaction strengths, the radiative decay rate of an 
emitter is directly proportional to the projected local density 
of states at the position ro of the emitter, i.e., 7 ra d(ro, A) oc 
|p(A)| 2 p(ro, A), where p(A) is the transition dipole moment, 
and A is the wavelength of the optical transition. Thus the 
QDs embedded in the samples are direct probes of p(ro) and 
from its spatial fluctuations the Co correlation is extracted. 
To that end, we measure single QD emission decay curves by 
time-resolved spectroscopy, as shown in Fig. |3^a), and per- 
form the ensemble average by repeating the measurements 
for QDs placed at different positions. The measured decay 
curves are biexponential revealing the fine structure of the 
lowest exciton stated It splits into bright and dark exciton 
states, which differ in their total angular momentum being 
lh or 2h, respectively. For this reason, a dark exciton cannot 
recombine through an optical dipole transition unless it be- 
comes bright after a spin-flip process^-'. From the biexponen- 
tial model we extract the spin-flip rate between the dark and 
bright excitons, the non-radiative, and the radiative decay 
rate from the bright exciton of different QDs, as explained in 
detail in Ref. [19]. According to our measurements, only the 
radiative decay rate is affected by disorder while the other 
decay rates remain unperturbed (not shown). 

Figure [3jb) displays the histogram of 7 ra d measured in an 
unpatterned substrate (bulk) and in photonic crystals with 
different amounts of disorder. All the decay curves are mea- 
sured in a narrow wavelength range (A = 910±3nm) in order 
to avoid any effects of the wavelength dependence of the tran- 
sition dipole momen1r= For no intentional disorder (S = 0%), 
we observe a strong suppression of the spontaneous emission 
rate compared to bulk due to the 2D photonic band gap. 
With increasing amount of disorder, the distribution of de- 
cay rates broadens tremendously, cf. Fig. [3jb), and largely 
enhanced rates are observed in the regime of Anderson lo- 
calization (8 > 6%). For S = 6% in particular, we observe a 
maximum rate of 7 ra d = (9.13 ± 0.02) ns _1 and a minimum 
rate 7 ra d = (0.08±0.06) ns _1 , see FigpTa), which corresponds 
to a Purcell enhancement as high as 114. To our knowledge 
this is the largest Purcell enhancement ever reported in any 
nanophotonic quantum electrodynamics experiment, and il- 
lustrates the great potential of disordered photonic materials 
for light-matter enhancement. 

Figure|3|c) compares the measured Co correlation function 
to the intensity fluctuations of the Anderson-localized modes 
as recorded from the far-field measurements shown in Fig. [2] 
after averaging over the same wavelength range as we mea- 
sured 7rad. A close similarity is observed, which is surprising 
since the intensity-fluctuation measurements probe the local- 



ized modes after excitation by many QDs, which is not a 
direct probe of p(i"o) as required for Co measurements. Hav- 
ing extracted both g and Co independently enables quanti- 
fying the nonuniversality of the latter. In the regime where 
scattering can be treated perturbatively Co is inversely pro- 
portional to y/p, and deviation from this scaling is a signature 
of nonuniversality^^. In Fig. j^d) we plot ^fg x Co, which 
clearly deviates from the universal scaling as a result of the 
breakdown of the perturbative scattering theory in the local- 
ized regime and the fact that Co depends on the near-field 
properties of the source. Consequently Co is sensitive to the 
near-field environment of the source, i.e., information about 
the correlation length of the local dielectric function could 
be extracted. This sensitivity may pave a route to subwave- 
length diffuse- wave microscopy of disordered opaque photonic 
medijP. Finally we have observed that the Co correlation 
function can be tuned by taking advantage of the underlying 
order of the photonic crystal lattice. In Fig. |3|b) we show 
for S = 6% the decay rate distributions for two disordered 
photonic crystals with the same lattice constant but differ- 
ent hole radii. By changing the hole radii we can tune the 
photonic band edge to be out of resonance with the quantum 
dot luminescence spectrum (corresponding to yellow data of 
Fig.|3|b)) thereby significantly reducing the strength of multi- 
ple scattering and hence Co- We obtain Co = 0.62(8), which is 
five times smaller than the value measured at the band edge, 
revealing that the Co correlations can be efficiently tuned by 
controlling the dispersion. 

We have presented the experimental observation of nonuni- 
versal correlations when embedding light emitters in a com- 
plex random medium. The ability to measure these corre- 
lations opens the possibility of probing locally the complex 
properties of random media. Such insight is essential for us- 
ing random media for lasing^or quantum optics experiments 
where potentially a single photon and a single quantum dot 
can become entangled by exploiting disorder for enhancing 
light-matter interaction*^. The nonuniversality of the Co cor- 
relations may be exploited for superresolution spectroscopy 
inside an opaque medium, since it gains insight into the near- 
field properties. Random media are currently also considered 
as a way to harvest light for solar energy applications by 
employing enhanced local absorptiorP^, and our work and 
method may be applied for engineering the efficient absorp- 
tion of light. 
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SUPPLEMENTARY INFORMATION 

Fabrication of the photonic crystal structures. 

The samples are GaAs membranes containing a single layer 
of self-assembled InAs quantum dots (QDs) in the center 
with density 80 |j.m~ 2 . An ordered triangular lattice of holes 
is etched in the structure creating a 2D photonic band gap 
that suppresses in-plane propagation of light. We investi- 
gate a set of samples with dimensions 7 um x 7 um, mem- 
brane height h — (150 ± 5) nm, a range of lattice constants 
220 nm < a < 280 nm, and hole radii 0.24a < r < 0.3a. 
The hole positions are randomly varied using a Gaussian ran- 
dom number generator function. The amount of disorder 5 
is characterized by the standard deviation of the hole posi- 
tion with respect to the lattice constant varying from S = 0% 
to 5 — 12%. The photonic band structure was calculated for 
8 = 0% using the software MIT PHOTONIC BANDS. The 
refractive index of GaAs, n = 3.55 — 3.44, is calculated^ for 
T = 10 K in the spectral range A = 850 nm — 1000 nm. 

Photoluminescence experiments. 

For optical characterization we use a confocal micro- 
photoluminescence setup to excite QDs within a diffraction- 
limited region and collect the emitted light (details can be 
found in Ref. |S2|). The samples are cooled to a temperature 
of T = 10 K in a helium flow cryostat and the sample posi- 
tion is controlled with stages with a resolution of 0.3 um. To 
efficiently excite Anderson-localized modes, we employ the 
inhomogeneously broadened spectrum from saturated QDs 
pumped at a high excitation power density of P = lkW/cm 2 . 
Emission spectra are collected in a wide wavelength range 
850 nm — 1000 nm with a spectrometer (0.15 nm resolution) 
equipped with a CCD camera while scanning the excitation 
and collection objective along the photonic crystals. The in- 
tensity probability distribution P(I) is measured by collecting 
the intensity I^, y at each spatial position (a;, y) with a spatial 
binning size of 0.23 um and at different wavelengths with a 
binning size of 1 nm. Subsequently, an average over the wave- 
length range 850 nm — 1000 nm is performed. 

The dimensionless conductance, g, is obtained by fitting 
the experimental intensity probability distribution with the 
full speckle intensity distribution obtained in Ref. [S3 in the 
regime of perturbative scattering and the absence of absorp- 
tion as: 



P(I) 



Ax 



(2) 



where Ko is a modified Bessel function of second kind and 
&con(x) is obtained by assuming plane- wave incidence as: 



$con(a:) = gln{Jl + - + J -) 



(3) 



Probing the radiative decay rate. 

The QD decay rate is measured by time-resolved photo- 
luminescence spectroscopy at an excitation power density of 
P = 24W/cm 2 , which is below the saturation level of a sin- 
gle QD. The QDs are excited with short optical pulses (2ps) 
at a 75 MHz repetition rate. Single QD lines are detected by 
spectral filtering and their decay curves are recorded with 
an avalanche photo diode. By accounting for the exciton fine 
structure, it is possible to extract the radiative decay rate and 
therefore to separate effects from nonradiative recombination 
and spin-flip processes. The decay curves are fitted with a bi- 
exponential model I(t) — Ate 'f t + A s e ,st , where subscripts 
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s and / denote slow and fast decay components, respectively. 
From 7/, j a and Af/A s we determine the spin-flip rate be- 
tween the dark and bright excitons, the non-radiative, and 
the radiative decay rate of the bright excitons by employing 
the method detailed in Ref. |S4| . We note that the proper ex- 
traction of the radiative rate, as opposed to merely studying 
the fitted rates 7/ or y s , is essential in order to probe the local 
density of optical states and thereby obtain the Co correlation 
function. 
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